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ABSTRACT 

Hypervelocity  flows  for  velocities  in  excess  of  1.4  km/sec  (Mach  5)  require 
very  high  stagnation  temperature  to  avoid  liquefaction.  The  arc  heater  wind  tunnel 
has  been  designed  to  provide  such  flows.  The  electric-arc  driven  wind  tunnel  can 
develop  stagnation  temperatures  up  to  13,0(X)''K  which  will  produce  hypervelocity 
flows  up  to  7  km/sec  (earth  orbital  speed).  The  nature  of  the  flow,  however,  is  such 
that  the  high  temperature  source  flow  may  cause  severe  gradients  at  the  nozzle  exit. 
In  order  to  perform  aero-thermodynamic  tests  the  characterization  of  the  flow  in  the 
test  section  is  required.  This  paper  experimentally  determines  the  stream  profiles  for 
an  arc  jet  wind  tunnel  conical  nozzle  direcUy  from  calorimetry  and  pitot  probe 
surveys. 
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I.    INTRODUCTION 

This  report  is  concerned  with  developing  a  procedure  by  which  flow  field 
characteristics  of  very  high  temperature  arcjet  flows  may  be  produced  directly  from 
stream  surveys  of  nozzles. 

Hypervelocity  flight  in  the  atmosphere  has  previously  been  concerned  with 
entry  into  the  atmosphere  and  descent  to  the  ground.  This  has  been  expanded 
during  the  last  decade  to  include  vehicles  which  will  use  the  upper  layers  of  the 
atmosphere  to  perform  hypervelocity  aerodynamic  maneuvers  and  vehicles  which 
will  fly  from  the  ground  into  orbit  using  air  breathing  powered  lift.  The  design  and 
testing  of  an  aeroassisted  orbital  transfer  vehicle  (AOTV)  and  the  National  Aerospace 
Plane  (NASP)  will  require  ground  based  testing  in  order  to  obtain  a  data  base  for 
this  flight  regime.  Such  data  will  be  used  to  validate  3-D,  reacting  real-gas,  CFD 
codes  that  will  be  used  in  the  design  of  advanced  vehicles  where  aero- 
thermodynamic  effects  are  important. 

The  60-MW  Shuttle  Interaction  Heating  Facility  (IHF)  was  designed  to  meet  the 
testing  requirements  for  development  of  the  thermal  protection  system  of  the  shuttle 
orbiter  vehicle.  Its  design  provided  for  heat  transfer  simulation  and  stagnation- 
region  flow  simulation  for  conditions  which  the  shuttle  orbiter  would  encounter.  If 
this  facihty  is  to  be  used  in  aerodynamic  testing,  the  arcjet  flow  characteristics  need  to 
be  known  in  order  to  determine  if  they  provide  optimum  simulation  of  the  actual 
flow  conditions. 

In  arcjet  operations,  extremely  high  temjjeratures  are  generated.  For  example, 
for  the  test  reported  here,  the  stagnation  temperatures  are  of  the  order  of  6000°K 
(10,800°R).    Probes  must  be  highly  cooled  in  order  to  survive  the  environment. 
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Conventional  techniques  using  thermocouples,  for  example,  to  determine  local 
stream  energy  simply  will  not  work  because  of  these  high  temperatures.  Stream 
tube  flow-swallowing  techniques  can  be  used  up  to  moderate  temperatures, 
however,  such  enthalpy  probes  require  sharp  leading  edges;  and  these  can  not  be 
cooled  in  the  arcjet  environment.  Moreover,  such  probes  with  blunt  cooled  leading 
edges  yield  an  uncertain  stream-tube  capture;  and  resolution  and  accuracy  both  have 
large  uncertainties.  To  circumvent  these  problems  for  arcjet  stream  surveys,  a  simple 
method  has  been  developed  that  uses  small  diameter,  cooled,  hemispherical-shaped 
probes  to  measure  local  heat  transfer  and  pitot  pressure.  It  is  shown  in  this  report 
that  these  two  measurements  yield  stream  profiles  to  a  high  degree  of  resolution. 
The  theoretical  method  derived  here  has  been  formulated  into  a  computer  code  for 
data  reduction  of  stream  surveys.  This  code  can  be  incorporated  into  a  real-time  data 
acquisition  system  and  will  provide  on-line  reduced  stream  profiles  during  the 
nozzle  survey  for  immediate  assessment. 

PROFILE  is  a  computer  code  developed  to  use  the  60- MW  IHF  stream  survey 
data  for  a  13-inch  conical  nozzle.  The  PROFILE  code  directly  produces  normalized 
velocity,  density,  mass  flux  and  enthalpy  profiles.  The  derivation  of  PROFILE  uses 
the  basic  conservation  equations  and  is  described  in  Chapter  IV.  PROFILE  uses  the 
heat  transfer  and  stagnation  pressure  data  as  inputs  and  outputs  normalized  values  of 
heat  transfer  rates,  stagnation  pressure,  velocity,  density,  mass  flux  and  enthalpy 
with  respect  to  centerline  values  of  the  flow.  Characteristics  of  flows  for  other  arcjet 
designs  and  other  nozzles  can  be  readily  obtained  using  heat  transfer  and  stagnation 
pressure  survey  data  in  conjunction  with  the  PROFILE  code  developed  here. 


II.     INTERACTION  HEATING  FACILITY  DESCRIPTION 

A  schematic  diagram  and  photograph  of  the  IHF  arcjet  facility  are  shown  in 
Figure  1  and  2,  respectively.  The  major  components  of  the  facility  include  a 
constricted-arc  air  heater,  a  150  NfW  power  supply,  two  interchangeable  hypersonic 
Mach  number  nozzles,  a  walk-in  test  chamber,  and  ancillary  subsystems  consisting 
of  a  steam-ejector  vacuum  system,  cooling  water  system,  and  automated  data 
acquisition  system.  The  high-power  capability  of  the  arc  heater  (rated  to  60  MW) 
provides  for  full  scale  subsystem  tests  in  high  enthalpy  streams  in  both  large  scale 
stagnation  flows  and  boundary  layer  flows  at  high  Reynolds  numbers. 

The  primary  component  of  the  facility  is  a  constricted-arc  heater  390  cm 
(153.3  in.)  long  and  8  cm  (3.15  in.)  in  diameter  shown  in  Figure  3.  The  rated 
airflow  rate  is  1.4  kg/s  (3.09  Ibs/s)  and  the  rated  pressure  is  12  atm.  Cooling  water 
is  supplied  at  7000  kPa  (1(XX)  psi),  which  provides  cooling  for  a  maximum  heating 
of  12  kW/cm-  to  the  constrictor  wall.  Eight  anode  and  cathode  electrodes,  each 
capable  of  conducting  1000  A,  provide  a  maximum  arc  current  capacity  of  8000  A. 
Thus,  the  capability  for  continuous  operation  at  high  current  and  high  enthalpy 
levels  is  provided. 

Direct  current  for  operating  the  arc  heater  is  provided  by  a  150  MW  D.C. 
power  supply.  Six  modules  make  up  the  power  supply,  each  consisting  of  a  three- 
phase,  full-wave,  phase  controlled  silicon  rectifier.  One  module  is  capable  of 
operating  up  to  a  level  of  5500  V  into  an  arc  load  at  current  levels  up  to  27(X)  A. 
The  six  modules  can  be  connected  in  any  compatible  series/parallel  arrangement  by 
means  of  remotely  operated  setup  switches.  The  duty  cycle  for  the  power  supply  is 


30  minutes  on,  30  minutes  off.  A  150  MW  rating  is  achieved  for  run  durations  of 
approximately  15  seconds. 

The  facility  has  two  interchangeable  nozzles:  1)  a  semielliptic,  Mach  number  5 
nozzle  with  an  open-jet,  flat-plate  test  section  measuring  80  X  80  cm  (31.5  X  31.5 
in.)  for  boundary-layer  surface  flow  tests  and  2)  a  conical,  nozzle  for  testing  large- 
scale  stagnation  flows  in  a  free  jet.  The  conical  nozzle  is  constructed  in  5  sections 
with  exit  diameters  varying  from  15.24  cm  (6  in.)  to  105  cm  (41.4  in.).  The  Mach 
number  varies  with  the  nozzle  area  ratio  from  3  to  7.5  depending  on  the  exit 
diameter  selected.  Both  nozzles  are  water  cooled  by  intemal  water  passages,  and  the 
throat  diameters  are  6.03  cm  (2.37  in.).  Both  nozzles  are  designed  to  flow  into  an 
evacuated,  walk-in  test  chamber  that  can  be  maintained  at  a  pressure  of  less  than 
0.004  atm  (3  mmHg),  a  pressure  that  varies  with  flow  rate  according  to  the 
pumping  curve  of  the  steam  ejector  system.  Tests  have  shown  that  supersonic  flow 
with  33  %  blockage  ratio  has  been  attained  in  the  conical  nozzle,  corresponding  to  a 
model  diameter  of  about  60  cm  (24  in.). 

An  on-line  automated  data  acquisition  and  computing  system  is  used  to 
record  arcjet  performance  and  test  data.  The  conversion  from  analog  to  digital 
signals  is  made  by  means  of  an  integrating  digital  voltmeter,  so  that  the  signals  can  be 
processed  by  an  on-line  computer.  Output  consists  of  a  real-time  display  of  five 
channels  that  are  monitored  during  a  run,  printed  and  reduced  data,  and  a 
continuous  time  history  data  plot.  [Ref .  1  ] 


III.     SURVEY  DESCRIPTION 

The  design  of  the  IHF  allows  for  the  control  of  two  variables:  1)  constrictor 
pressure  and  2)  arc  current.  The  test  conditions  are  uniquely  determined  by  these 
two  control  settings  only.  For  this  survey,  the  test  chamber  pressure  was 
maintained  at  a  low  pressure  to  allow  the  measurements  to  be  made  in  a  free  jet.  No 
attempt  was  made  to  control  the  expansion  ratio  into  the  free-jet  chamber.  In  general, 
the  flow  was  as  an  under-expanded  free  jet.  This  type  of  operation  causes  the  free  jet 
to  spread  and  to  operate  at  a  slightly  higher  Mach  number  relative  to  the  nozzle  exit 
plane.The  survey  measurements  close  to  the  nozzle  exit  plane,  then  are  relatively 
uninfluenced  by  the  free-jet  expansion  process.  Most  of  the  discussion  pertains  to 
the  7.62-cm  (3-inch)  station  to  avoid  the  complexity  of  the  additional  expansion. 
Constrictor  pressure  was  varied  for  several  values  of  current  in  order  to  map  the 
operating  envelope  for  stagnation  flows  using  the  conical  nozzle. 

The  probes  used  to  obtain  data  consisted  of  a  hemispherical  Garden  type 
calorimeter  and  a  pitot  probe  shown  in  Figures  5  and  6.  The  probes  were  1.58  cm 
(0.625  in)  in  diameter  and  were  mounted  20.32  cm  (8.0  in)  apart  horizontally.  The 
probes  were  mounted  to  a  transverse  mechanism  that  was  positioned  to  pass 
horizontally  through  the  flow.  Sweep  measurements  were  taken  at  7.62  cm  (3.0 
in),  19.5  cm  (7.7  in),  34.8  cm  (13.7  in)  and  52.6  cm  (20.7  in)  downstream  of  the 
nozzle  exit. 

Figure  7  shows  probes,  probe  supports  and  their  relationship  to  the  nozzle 
during  a  test  using  a  15.24  cm  (6.0  in.)  exit  diameter  nozzle.  Figure  8  shows 
probes  entering  the  flow.   Figure  9  shows  probes  approximately  on  centerline  of 


nozzle.  All  photographs  were  taken  using  only  radiated  light  from  gas  behindthe 
normal  shocks  on  probes  and  support  mechanisms. 

Data  consisted  of  heat  transfer  rate  and  stagnation  pressure  measurements  taken 
at  0.212  cm  (0.083  in.)  intervals  across  the  flow.  Typical  profiles  for  heat  transfer 
rate  and  stagnation  pressure  for  each  current  setting  are  shown  in  Figure  10. 


IV.    DERIVATION  OF  DIRECT  METHOD  FOR  PROFILES 

Determining  the  stream  profiles  for  an  arcjet  wind  tunnel  nozzle  requires  that  the 
heat  transfer  rate  and  the  stagnation  pressure  measurements  behind  a  nonnal  shock 
make  up  the  only  parameters  on  which  the  velocity  characteristics  depend.  Briefly, 
this  method  is  based  on  the  fact  that  stagnation-pwint  heat  transfer  is  proportional  to 
the  product  pV^  (approximately)  while  the  pitot  pressure  is  proportional  to  pV^. 
When  these  two  quantities  are  known,  a  simultaneous  solution  gives  density  and 
velocity  directly  from  the  measurements.  Once  the  velocity  profile  has  been 
determined  the  profiles  for  enthalpy,  mass  flux  and  density  can  be  readily  found. 

In  the  following  derivation,  free  stream  flow  parameters  are  those  found  in  front 
of  the  normal  shock  and  will  be  designated  by  subscript  one  while  subscript  two 
will  denote  parameters  behind  the  normal  shock.  The  starting  point  is  the 
momentum  and  continuity  equations  in  the  following  form: 

Pl^l=P2^2  (2) 

Rearranging  the  momentum  equation  and  substituting  into  it  the  continuity  equation 
gives: 
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Solving  for  p2 
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Noting  that  M2  is  much  less  than  1 ,  that  M2  =  (p(  Y,  M)  and  recalling  that  from  state  2 
to  stagnation  state  2  is  along  an  isentrope  then: 


^=''2(l^(^>2) 


t/(y-i) 


(5) 


Normalizing  equation  (5)  with  centerline  values  results  in 
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Substituting  equation  (6)  into  a  nomialized  equation  (4)  gives: 
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At  this  point  the  normalized  stagnation  pressure  behind  the  shock  is  related  to 
the  velocity  in  front  of  the  shock.  The  next  step  introduces  the  closed  form  of  the 
Fay-Riddell  [Ref  2]  equation  derived  in  Appendix  A.,  which  is  a  correlation  for  the 
heat  transfer  rate  to  a  hemisphere;  namely, 


*  =  «^^1^7^(*--'V)(Btu/sec-ft2)        (8) 

where  ht  is  the  total  enthalpy  of  the  flow  and  hw  is  the  enthalpy  of  the  wall.  The 
units  are:  h,  Btu/lbm;  pt2,  atm;  and  R,  ft.  Normalizing  with  centerline  values  gives: 
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Using  the  following  form  of  the  energy  equation 
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and  equation  (7)  in  equation  (9)  while  neglecting  the  enthalpy  of  the  wall  compared 
to  the  total  enthalpy  of  the  flow  yields: 
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Since  the  velocity  in  front  of  the  shock  is  supersonic  the  last  term  drops  out. 
Changing  notation  to  reflect  only  centerUne  and  free  stream  gives: 
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Squaring  equation  (13)  and  dividing  by  equation  (7)  removes  die  density  term  and 
results  in  an  equation  which  relates  free  stream  velocity,  stagnation  pressure  and  heat 
flux. 
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This  equation  equates  the  normalized  velocity  of  the  free  stream  to  the  two  measured 
quantities  behind  the  normal  shock. 

The  assumption  made  to  derive  this  result  are  simply  that  the  flow  is 
axisymmetric  and  the  enthalpy  of  the  wall  is  small  compared  with  the  total  enthalpy 
of  the  free  flow.  The  wall  enthalpy  is  measured  directiy  from  cooling  water  and  is 
found  to  be  approximately  0.3  MJ/kg  (134  Btu/lbm)  or  less  than  5  %  of  total 
enthalpy. 

Once  the  normalized  velocity  has  been  found  it  is  used  in  equation  (7)  to  find  the 
value  for  the  normalized  density. 


Pt 


'CL 


''CL^'^CL 


(15) 


Normalized  mass  flux  is  obtained  from  the  product  of  normalized  density  and 
velocity. 


mf    ^_P 
mL.        P 


'CL 


CL   kYcL 


(16) 


The  normalized  enthalpy  is  taken  directly  from  the  energy  equation  where 
gravitational  effects  are  neglected  and  the  static  enthalpy  is  assumed  to  h>e  much  less 
than  total  enthalpy,  which  results  in  the  following  equation. 
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*CL      VclJ  (17) 

Equations  (14),  (15),  (16),  and  (17)  were  used  in  the  code  PROFILE  to  calculate  the 
characteristic  profiles  of  the  arcjet  high  enthalpy  flows  using  survey  data. 


11 


V.    RESULTS 

For  high-powered,  arc-heated,  hypersonic  facilities  the  testing  for  aerodynamic 
or  thermodynamic  (heat  transfer)  studies  requires  a  central  core  flow  with  properties 
well  characterized.  Ideally,  of  course,  constant  core  flow  properties  are  desirable. 
The  present  study  was  undertaken  to  characterize  the  flow  properties  of  the  60-MW 
arc  heated  wind  tunnel  at  NASA- Ames  Research  Center.  Such  flow  properties  are 
essential  to  the  engineer  for  interpretation  of  the  results  for  aero-thermodynamic 
tests. 

Arc  heaters  operate  with  a  large  enthalpy  gradients  near  the  walls.  These 
gradients  cause  the  flow  properties  in  the  central  core  to  vary;  and  are  an  important 
concern  in  simulation  in  these  ground-based  test  facilities.  The  degree  of  this 
variation  has  not  been  reported  up  to  the  present  time  for  the  arcjet  facility  at  NASA- 
Ames  Research  Center.  Enthalpy  gradients  are  inherent  in  all  types  of  arc  heaters 
that  operate  with  cooled  walls.  For  example,  the  temperature  ratio  of  the  wall  for  a 
typical  arc  heater  varies  from  2  -  5  %  when  the  enthalpy  level  is  in  the  typical  vehicle 
reentry  range  of  18.6  -  41.87  MJ/kg  (8000-18000  Btu/lbm)  which  in  terms  of 
velocity  equates  to  20,(XX)-30,CKX)  ft/sec.  The  constricted  arc  heater  was  developed 
to  minimize  this  gradient  by  the  current  density  and  length  to  cause  the  enthalpy 
profile  to  assume  a  flatter  shape  than  a  free-burning  or  vortex-stabilized  arc  heater. 
The  theory  of  Stine- Watson  [Ref.  3]  shows  that  for  a  constricted  arc  heater  the 
enthalpy  profile  assumes  a  Bessel  function  form  as  opposed  to  a  linear  profile  for  a 
free-burning  or  vortex-stabilized  design  (Huels  arc)  illustrated  in  Figure  1 1 . 

Results  from  the  profile  measurements  obtained  here  can  be  used  to  answer  the 
question  as  to  whether  the  property  variables  (enthalpy,  velocity,  density  and 
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energy  flux)  have  the  same  profile  or  shape  at  the  nozzle  exit  as  they  do  for  the 
constrictor.  Scant  measurements  exist  for  the  constrictor.  Most  of  the  values  for  the 
constrictor  arc  heaters  are  predicted  by  the  Stine- Watson  theory.  Elaboration  of  this 
theory  is  represented  by  a  computer  code  called  ARCFLO  n  [Ref.  4].  This  code 
was  developed  to  predict  the  performance  of  high-pressure,  high-enthalpy 
constricted-arc  heaters.  The  conservation  equations  (mass,  axial  momentum,  and 
energy)  are  written  for  a  steady-state  and  radiatively  participating  gas  mixture  with 
electric  current  conduction.  The  governing  equations  employ  the  boundary-layer 
approximation  and,  hence,  the  pressure  in  the  wall  normal  direction  is  assumed 
constant.  Auxiliary  relations,  such  as  the  equation  of  state,  global  current  continuity, 
and  global  mass  continuity,  are  used  to  relate  the  gas  density  with  the  pressure,  mass 
flow  rate,  and  mixture  properties. 

The  profiles  of  density,  velocity,  enthalpy  and  energy  flux  computed  by 
ARCFLO  n  in  the  constrictor  are  compared  with  experimental  survey  data  obtained 
at  the  nozzle  exit  in  Figure  12.  The  difference  seen  in  the  density  profiles,  Figure 
12a,  can  be  accounted  for  by  the  increase  in  density  near  the  wall  of  the  constrictor  as 
opposed  to  the  decrease  in  density  in  an  expanding  free  jet.  The  velocity  profiles. 
Figure  12b,  are  very  similar  with  the  experimental  data  showing  a  slightly  flatter 
profile.  The  enthalpy  profiles,  Figure  12c,  shows  for  ARCFLO  II  a  central  core 
region  where  the  enthalpy  drops  moderately  from  the  centerline  to  r/R  =  0.4,  a 
linearly  decreasing  region  to  r/R  =  0.95,  and  then  a  steep  drop  near  the  wall  where 
the  conduction  losses  dominate.  The  ratio  of  centerline  enthalpy  to  mass  average 
enthalpy  is  1.6.  This  low  ratio  value  is  typical  for  wall-stabilized  arcs.  The  survey 
data  exhibits  a  much  flatter  profile  than  ARCFLO  II  with  a  ratio  of  centerline 
enthalpy  to  mass  average  enthalpy  of  1.15. 
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The  product  of  the  above  three  profiles  is  shown  in  Figure  12d  and  is  the 
energy  flux  of  the  airstream.  Note  that  this  parameter  is  proportional  to  the 
convective  heat-transfer  rate  to  a  body  within  the  stream  and  that  it  is  a  measure  of  the 
usable  stream  area.  The  flat  profile  for  ARCFLO  n  implies  that  stagnation  point 
heating  to  a  test  body  would  be  constant  to  within  10  %  of  the  value  obtained  on  the 
centerline  flow  for  r/R  up  to  0.7.  The  survey  data  profile  shows  a  more 
pronounced  curve  which  can  be  attributed  to  the  decrease  in  density  in  the  free  jet.  If 
test  chamber  static  pressure  is  maintained  at  the  nozzle  exit  static  pressure  then  the 
density  profile  would  become  flatter  and  thereby  improve  the  energy  flux  profile. 

Figures  13  -  28  are  plots  of  PROFILE  code  output  data  at  specified  survey  test 
points.  All  test  point  plots  are  normalized  with  respect  to  the  corresponding  values 
obtained  at  the  3-inch  downstream  positions. 
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VI.     CONCLUSIONS  AND  RECOMMENDATIONS 

The  present  study  of  survey  data  for  the  13 -inch  IHF  nozzle  showed: 
That  there  is  good  agreement  between  the  ARCFLO II  code 
prediction  for  velocity  profile  and  experimental  data; 
Enthalpy  profile  of  experimental  data  is  much  flatter  than  ARCFLO  n 
theory  predicts,  which  can  be  explained  by  the  codes  inability  to 
account  for  energy  transfer  toward  constrictor  edges. 

It  is  recommended  that: 

Future  surveys  control  test  chamber  conditions  in  order  to  determine 

if  improvements  in  the  profile  characteristics  result; 

Sting  mounted  heat  transfer  and  stagnation  pressure  probes  be 

developed  in  order  to  reduce  probe  mount  interference  with 

boundary  layer  of  the  flow; 

Survey  data  be  organized  into  accessible  files  so  that  characteristic 

maps  for  individual  nozzles  and  complete  nozzles  can  be  utilized  by 

researchers  and  contractors. 
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APPENDIX  A 

DERIVATION  OF  THE  CLOSED  FORM  OF  THE  FAY-RIDDELL  EQUATION 

by  WARREN  WINOVICH 

1.  The  Nusselt  number  is  defined  by: 


H  X         / 
Nu  =  -r —        (w  =  wall  value) 


2.  The  heating  rate  is  given  by: 


where 


Cp  -        h,-h^ 

C    "=■  c      =  =  1 

Cp  ^p       T   -T 


(h.-h^) 


where  subscript  "e"  indicates  external  flow  conditions 
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(Al) 


^'-  ^  (A2) 


(by  definition  of  Nusselt  number;  wall  temperature  used  as  reference) 
3.  Solving  for  heating  rate  in  terms  of  Nusselt  number. 


^'-  ^  (A3a) 


(A3b) 


'-V^V^^/P^-P^:^  (A3C) 


4.  For  the  laminar  boundary  layer  near  the  stagnation  point: 


(A4a) 


therefore 


(A4b) 


5.  Evaluation  of  velocity  gradient  at  the  stagnation  point  for  a  hemisphere: 


/?,  + -r- -p^ «,^   =  p^         (momentum  equation) 


2gc  P, 


Pe 

P, 

1 


(A5a) 

p  /  X  \ 

-p-^  =  1 -sin^f -^j         (Newtonian    flow  over  hemisphere) 


''■-M)Ht) 


(A5b) 


dx  J      /?  V    ^'yP.  J        V/?/  (A5c) 
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dx  J        R 


'2g 


r  P.  \ 


vP',y 


1  P ,  =  P,      for     X  =  0 


(A5d) 


6.  The  heating  rate  equation  becomes: 


Nu    I —  1    /TT^ 


{K-K) 


Pr,  c 


(A6) 


7.  The  Fay  and  Riddell  result  for  the  heating  transfer  coefficient  is: 


«.4 


^«     =0.763(/>rn#4^1    c 


Pw^A. 


(equation  63  and  45  Fay  and  Riddell) 

8.  Substituting  equation  (A7)  into  equation(A6);  using 


(A7) 


f  Pt  \ 


V^'W 


^R  T,Z,         and     R=5a5 


gives: 


q=   0.763(Pr)''V2gc(53.3)Vp7]i; 


Pwji.. 


e.4 


-^fTzT, 


Vr 


(A8a) 


4^^        0.763    ,/i — rrr^  nrTr(^j:}!^\\rrT' 


(A8b) 
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9.  Using: 


p,  =  39.65j,  J 


•  .<2S 


H.  =  184xl0-'(j^)       K. 


(A9a) 


(A9b) 


where 


w,   Keitct 


^  fornuU 

Equation  (A9b)  is  appropriate  for  temperatures  of  6000*K  (lOOOO'R)  where  arcjet 
operate  and  for  reentry  and  reentry  stagnation  conditions  (Vgntry  >  5000  m/sec  ~ 
15000  ft/sec) 


^  3 9.65 X  184  X  10"'  PJ^ 
1000  ^'^' 


^  ••*28 

^         =9.60x  10  *     •         K 


P.  '  T,Z,       *  (A9c) 


•  .(28 


p.^i.         T,Z,  T^  K 


^»  (A9d) 

Note  that  Kw  =  1 .00  for  the  expression  chosen  for  viscosity. 

10.  The  final  form  for  the  expression  for  heating  rate  becomes: 


-rr7-\/2gc(5  3.3)  3.10X  10  ^ 


VTsih^-K)-  (p,)-^    ^'^     -  '    '  ^T^z,) 


o.s 


•  .1  •.•i2S  t.l 

^''rlzT    IT"        Ik" 


(AlOa) 
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Simplifying  for  Pr  =  0.71  the  equation  becomes: 


=  {constant  )  - 


y/p~i{h,-  hj) 


(T^Z^) 


•.•372 


^  =  {constant}    ^  '    '..,,,  1^1       K. 


■y/p7{h,-  h^) 


Noting  that  Zw  =  1.00 

For  Pr  =  0.71  the  constant  in  brackets  above  becomes: 


•.•372    I    2, 


{constant  }  = -^^^^  ^2g,(53.3)  3.10x  10"'         =  0.02225 

(Pr) 

and  the  expression  for  heating  rate  is: 


'V^ 


-s/pTih,-  h^) 


(  •••372  •.•<38  ^ 


(AlOb) 


The  term  within  the  braces  is  only  a  function  of  the  stagnation  state.  For  usual  wind 
tunnel  tests  the  stagnation  state  is  within  the  hmits: 

0.01  <  Ps  <   1  atmosphere  130  <  hs  <  20,000  Btu/lbm 

For  this  range  of  pressure  and  enthalpy  the  term  in  the  braces  is  essentially  a 
constant: 


i 


K\ 


0.0372      ,  ^•.•(38 

,JT,Z   \  {T,Z,) 


Ta. 


•  .251      K. 


=  1.873±  0.042 
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Therefore  the  hemisphere  formula  for  heating  rate  is: 


or: 


iVr 


-y/plih.-   K) 


=  {0.02225}{1.873}     =     0.0417 


q  =0.0417    I ^{h,-  /I.) 


where 


q  =  Btu/s    ft' 

/?,  =  atm     -^   p^  (measured) 

R  =  feet 

h  =   Btu/lb„ 

ID 

Symbols 

Cp 

Specific  heat 

h 

Eniahlpy 

k 

Theniial  Conductivity 

Pr 

Prandtl  Number 

R 

Radius 

T 

Temperature 

Z 

Compressibilty  Factor 

P 

Density 

^ 

Viscosity 

(AlOc) 


Subscripts 

e        External  Conditions 
s        Stream  Conditions 
t         Stagnation  Conditions 
w       Wall  Conditions 
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APPENDIX  B 

UNCERTAINTY  ANALYSIS  FOR  DERIVATIONS 

Uncertainty  analysis  of  the  heat  transfer  and  stagnation  pressure  probes  is  based 
on  manufacture's  claims  for  uncertainties.  In  the  case  of  the  heat  transfer  probe  the 
manufacturer  list's  accuracy  of  readings  between  3  &  4  %  of  the  value  recorded. 
For  the  following  analysis  the  value  of  5  %  is  used.  The  accuracy  of  the  stagnation 
pressure  probe  has  been  found  by  NASA  personnel  to  be  on  the  order  of  1/4  %  of 
the  recorded  value.  Using  these  values  the  uncertainty  of  the  normalized  velocity  is 
shown  below. 

Starting  with  equation  (14): 

V 


CL 


^  <1  ^  f  Pt  \ 


K^CLJ 


V        CLJ 


and  taking  the  logarithmic  differential  gives: 


Integrating 


In 


d  V 


CL 


CL 


=  |ln 


f   d    \ 


V^cl) 


"4  '"IP 


^d^_  d_^^ 

9  ^CL   ) 


CL 


1   {^_1ECL 

^yp      PcL 


V^    =    y  [(e  r)4  -  (e  r)9'cJ     -  4-  [(e  r)P  -{^  r)PczJ 
CL  ^ 

where  (cr)  is  uncertainty. 

In  the  worst  case  where  centerline  values  are  opposite  in  sign  of  the  original 
values  results  in 
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dV 


CL 


y  [(0.05) -(-0.05)]    --i  [(0.0025)  -(-0.0025)] 


4r-   =  4-  [OiO]   -4-  [0.005] 


CL 


SO  that 


d  V 


CL 


±  5  % 


Applying  the  velocity  uncertainty  to  the  other  jwoperty  equations  (15),  (16)  and 
(17)  gives: 


CL 


CL 


dp 


Ap 


d  V 

^cLn 


-1 


±  10% 


dm/ 


CL 


dp  V d  V 
P 


CL 


±  15  % 
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APPENDIX  C 

COMPUTER  PROGRAM  FOR  DETERMINING  PROFILES 

FILEi  PR0F2     FORTRAN   Al 

PROGRAM  PROFILE 

C  THIS  PROGRAM  IS  DESIGNED  TO  CALCULATE  FLOW  PROFILES  FOR 

C  NASA  AMES  IHF  ARCJET  HIND  TUNNELS.   VELOCITY,  DENSITY,  MASS  FLUX 

C  AND  ENTHALPY  PROFILES  NORMALIZED  TO  CENTERLINE  CONDITIONS 

C  ARE  DEVELOPED. 

DIMENSION  QDOT(IOOO),  QDOTN(IOOO),  PTOT(IOOO),  VNOR(IOOO), 
S  RHO(IOOO),  MFLUX(IOOO),  HTOT(IOOO),  DIST(IOOO), 

S  QSQR(IOOO),  PTOTN(IOOO) 

REAL  VNOR,RHO,MFLUX,HTOT,QSQR,PTOTN,QDOTN 

INTEGER  S,I,J,K,L,IMAX,IQD0T,IPT0T,VAL,I1,I2 

REAL  PMAX,  QMAX,VMIN 

CHARACTER  K36  G 

CHARACTER  XI  E,Y,N 

CHARACTER  X6  A,B,C 

CHARACTER  ^^5    F 

C       XXKXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXK)fXXXXXXXXXXXXXXXXXXXXXX 

C      THIS  SECTION  COPIES  TEST  CONDITIONS  INTO  OUTPUT  FILE 
C     xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
READ(5,1)F 

1  F0RMAT(A^5) 
WRITE(7,1)F 
READ(5,2)G,C 

2  FORMAT(A36,A6) 

WRITE(7,x)G,C  1 

DO  10  1=1,3 
READ(5,1)F 
HRITE(x,l)F 
WRITE(7,1)F 
10     CONTINUE 

C     xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
C      THIS  NEXT  SECTION  READS  TEST  DATA  FOR  QDOT  AND  PTOT. 
C      FIRST:  NORMALIZES  VALUES  OF  QDOT  AND  PTOT 
C      SECOND:  LOCATES  MAXIMUM  VALUE  TO  USE  AS  CENTERLINE 
C      THIRDi  CALCULATES  VALUES  FOR  VELOCITY,  DENSITY,  MASS  FLUX  AND 
C  ENTHALPY  PROFILES  NORMALIZED  TO  CENTERLINE 

C 

C      NOTEt  QDOT  AND  PTOT  DATA  SHOULD  CONTAIN  ONLY  ONE  ZERO  AS  LAST 
C  VALUE.   THIS  IS  CRITERIA  USED  TO  TERMINATE  READING  OF  DATA. 

C 

C     xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
IMAX=0 
DO  20  1=1,1000 

READ(5,15)QD0T(I),  PTOT(I),  DIST(I) 
IF(QDOT(I).EQ.0)GOTO  21 
IMAX=IMAX+1 

20  COflTINUE 

21  WRITE(X,x)«NUMBER   OF    DATA   POINTS    READ    =MMAX 
ENDFILE(5) 
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C      1  NORMALIZATION  OF  QDOT  AMD  PTOT 

C       XXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXKXXXXXXXXKXXXXXKXXKX 
QMAX  =  0 
PMAX  =  0 
DO  200  I  =  1,IMAX 

QMAX  =  MAX(QMAX,QDOT(I)) 
PMAX  =  MAX(PMAX,PTOT(I)) 
200   CONTINUE 

C       XXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXX 

C     2  DETERMINING  LOCATION  OF  MAX  VALUES.   THESE  WILL  BE  USED  TO 

C       ALIGN  PROBE  VALUES  TO  CENTERLINE. 

C     xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DO  201  I=1,IMAX 
IQDOT=I 

VAL=DIMCQMAX,QDOT(I)) 
IF(VAL.EQ.Q)GOTO  202 

201  CONTINUE 

202  DO  203  I=1,IMAX 

IPTOT=I 

VAL1=DIM(PMAX,PT0T(I)) 
IF(VAL1.EQ.0)G0T0  205 

203  CONTINUE 

205    WRITE(X,X)»I  AT  QMAX  =MQDOT,'   I  AT  PMAX  =',IPTOT 
C     xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
C      3  CALCULATION  OF  PROFILES 

C     xxxxxxxxx^xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
Il=IQDOT+179 
I2=IPT0T+179 

HRITE(X,X)'I1  =M1,  •     12  =M2 
K  =  0 

DO  300  I=IQD0T,I1 
K  =  K  +  1 

QDOTN(K)=ABS(QDOT(I)/QMAX) 
QSQR(K)=QD0TN(K)xx2 

300  DIST(K)=DIST(I)-DIST(IQDOT) 
K  =  0 

DO  301  I=IPT0T,I2 
K  =  K  +  1 

301  PTOTN(K)=ABS(PTOT(I)/PMAX) 
DO  302  1=1, K 

VNOR(I)  =  SQRT(SQRT(QSQR(I)/PTOTN(I))) 
RHO(I)  =  PT0TNC1)/(1/VN0R(I))XX2 
MFLUX(I)  =  RH0(I)XVN0R(I) 

302  HTOT(I)  =  VMQR(I)XX2 
GOTO  <400 

^00    HRITE(7,X)«QMAX  =  ',QMAX 
WRITE(6,X)'QMAX  =  • , QMAX 
HRITE(7,X)«PMAX  =',PMAX 
WRITE(6,X)'PMAX  =',PMAX 
IJRITE(7,x) •============== ========================================= 

8================! 

HRITE(7,X)'    HEATING     STAG     VELOCITY    DENSITY    MFLUX 
SENTHPY     DIST' 
WRITE(7,X)'      RATE      PRESSURE* 
WRITE( 7, X) •======================================================= 

HRITE(7,6  00)(QDOTN(I),PTOTN(I),VNOR(I),RHOCI),MFLUX(I),HTOT(I),DIS 
8T(I),I  =  1,K,'4) 
600   F0RHAT(6(5X(F5.3)),5X,F6.3) 
STOP 
END 
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Figure  8      Probes  entering  6-incli  nozzle  free  jet  flow  during  survey. 
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(a)  Heat  transfer  rate  verses  radius  at  3-inch 
position  from  nozzle  exit 
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(b)  Stagnation  Pressure  verses  radius  at  3-inch 
position  from  nozzle  exit 

Figure  10    Typical  profiles  of  the  heat  transfer  rate 
and  stagnation  pressure,  13-inch  nozzle  (A/A*  =  30) 
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Figure    11     Comparison  of  enthalpy  profiles  for  8-cm 
constricted  arc. 
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(a)  Density  Profile 


(b)  Velocity  Profile 
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(c)  Enthalpy  Profile 


(d)  Energy  Flux  Profile 


Figure    12     Comparison  of  ARCFLO  II  and  survey  data  profiles  of 
8-cni  constricted  arc. 
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